Multiple current reversals in forced inhomogeneous ratchets 
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j^JTJ' Transport properties of overdamped Brownian paricles in a rocked thermal ratchet with space 

^ | dependent friction coefficient is studied. By tuning the parameters, the direction of current exhibit 

multiple reversals, both as a function of the thermal noise strength as well as the amplitude of 
rocking force. Current reversals also occur under deterministic conditions and exhibits intriguing 
structure. All these results arise due to mutual interplay between potential asymmetry, noise, driving 
frequency and inhomogeneous friction. 
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Fluctuation induced transport in ratchet systems has been an active field of research over the last decade. In these 
systems in the absence of any net microscopic forces, asymmetric potential can be used to induce a unidirectional 
particle flow when subjected to an external nonthermal fluctuations. These studies have been motivated in part by 
the attempt to understand the mechanism of movement of protein motors in biological systems jl]]. To this effect 
several physical models have been proposed under the name of rocking ratchets , flashing ratchets |||| , diffusion 
*^ [ ratchets correlation ratchets ||, etc. In all these studies the potential is taken to be asymmetric in space. It 
has also been shown that one can obtain unidirectional current in the presence of spatially symmetric potentials. 
For these nonequilibrium systems external random force should be time asymmetric |l or the presence of space 
dependent mobility is required ]I^ , [io| , pT| , ^3| -[l^| . By suitably tuning the system parameters such as temperature, 
friction coefficient, mass, etc, one can even change the direction of the current. Indeed, the study of current reversal 
t-H , phenomena has given rise to a research activity on its own. The motivation being the possibility of new particle 
separation devices superior to existing methods such as electrophoretic method for particles of micrometer scale [fl7|| . 

Bartussek et. al. |Ij showed the occurrence of current reversal in rocked thermal ratchet with both amplitude of 
rocking force as well as the temperature of thermal bath. They attributed this current reversal to the "mutual interplay 
between noise and finite-frequency driving" . Multiple current reversals have also been shown in the deterministic 
limit of these ratchets when the inertial term is taken into account Ppp[|. However these multiple current reversals 
in inertial ratchets are not robust in the presence of noise. Beside rocking ratchets current reversals have also been 
observed in flashing ratchets (2|j2^,^|. We in this work report that that multiple reversals can be achieved even 
• in rocked overdamped ratchet in the presence of space dependent mobility, as a function of the noise strength and 
amplitude of the rocking force. More over these systems show current reversals when rocked adiabatically. In the 
Q deterministic overdamped case we get current reversal as a function of the amplitude of rocking force. Most of our 
results are attributed to the presence of space dependent mobility. We have studied the same system as due to 
Bartussek et al [I), except the presence of space dependent friction term. In earlier work the spatial asymmetry 
of the potential is responsible for unidirectional currents and their reversal as function of frequency. As mentioned 
^ , previouly one can get unidirectional currents in the presence of symmetric potentials, but for this space dependent 
friction is required p2| , ^0 ll,r^-|l6[. In these systems transport direction is given by the phase shift between mobility 
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and periodic potential. This phase shift induces spatial symmetry breaking as required for the directed motion. 
Appropriately choosing the phase shift leads to current reversal. We would like to emphasize that space dependent 
friction does not alter the equilibrium property of the system, however when the system is driven out of equilibrium 
not trivial dynamical effects arise due to space dependent friction p^JI^ , p^ | . Naturally in our present system we expect 
additional effects arising due to combination of spatial asymmetry and position dependence of friction coefficient. It is 
to be noted that systems with space dependent friction are not uncommon. Brownian motion in confined geometries 
show space dependent friction |2^] . Particles diffusing close to surface have space dependent friction coefficient [^3],|4| . 
It is believed that molecular motor proteins move close along the periodic structure of microtubules and will therefore 
experience a position dependent mobility |l5| ] . Frictional inhomogeneities are common in super lattice structures and 
Josephson junctions |^5| also. 

We consider an overdamped Brownian particle moving in a asymmetric potential V(x) with space dependent friction 
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coefficient r/(x) under the influence of external force field F(t) at temperature T. Throughout our analysis we take 
the ratchet potential V(x) = — ^(sin(27ra;) 4- j sin(47ra;)). Here \i is the asymmetry parameter with values taken in 
the range < [i < 1, friction coefficient rj(x) = r/o(l — Asin(27rx + 0)), |A| < f. <fi determines the relative phase shift 
between friction coefficient and potential . The forcing term is taken to be F(t) = Asin(wi + 0), (w = ^r, where r is 
period of force) . Without any loss of generality 9 is taken to be zero. The correct Langevin equation for this system 
in the overdamped limit is given by fl2|,^6|,p7| 



(V'(x) - F(t)) rfjx) k B T 

x = W) kBT WW + \W) m (1) 

where £(£) is a zero mean thermal Gaussian noise with correlation < >= 28{t — t') The equation of motion 

is equivalently described by the Fokker Planck Equation ( FPE) 

d -^r 1 = tt^t^I- + W - F ^ p M- (2) 

at OX T]{x) OX 

where P(x, t) is the probability density at position x at time t. Equation (||), in the form of a continuity equation 

dP{x,t) dJ(x,t) 



dt dx 

where 
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J(x,t) = -^-{(V>(x) F(t)) + k B T^-]P{x,t), (4) 
r](x) ox 

is the probability current. Since the potential and the driving force have spatial and temporal periodicity respectively, 
therefore J(x, t) = J(x + 1, t + t), pqfl. The average current j in the system is given by 
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lim - / dt J(x,t)dx. (5) 



It should be noted that for symmetric potential and A = 0, j = 0. Rectification of current is possible if the potential 
is either asymmetric or A ^ with (\> ^ 0, 7r fujfl . j is independent of the initial phase 9 of the driving force. 

We explore various parameter regimes of the problem extensively by solving the FPE, Eq. (|^) numerically with 
finite difference method. Throughout we have set current j and all other physical quantities such as T, A, w in 
dimcnsionless form . 

In the Fig. [j], the average current j is plotted vs T for various values of w. Here A = 0.1, (j) = 0.2tt, A = 0.5 and 
H=l. Unlike Bartussek et al. j|, where they show that current reversal is not possible under adiabatic conditions, 
however current reversal even in adiabatic condition can be obtained in the present case for some values of <fi []l6f . For 
moderately high frequency w = 4.0 the current reverses its sign twice at T = 0.1 and T = 0.22 and asymptotically goes 
to zero for higher values of T as shown in the figure. This phenomena of twice current reversal with temperature T is 
the foremost feature of our system, previously unseen in any overdamped system. The inset shows the zero contour 
plot of the current j versus T and <fi for three values of w. Crossing this zero contour line implies current reversal. It 
can be seen that twice current reversal occurs only for a very narrow range of 4>. For frequencies higher than certain 
critical frequency w c (4>) current flows in only one direction for all temperatures, i.e no current reversal occurs. The 
plot uu c vs 4> is shown in Fig. |2|. The adiabatic curve w < 1 is not shown in the figure as it goes much beyond the scale 
of the graph. However it has a similar qualitative shape as for w = 3.0 curve. As mentioned before, currents are due to 
the combined effect of phase shift (f> coming from space dependent friction and asymmetry parameter [x. In the regime 
4> = 0.2-7T and A = 0.5 and in the absence of asymmetry current flows in the negative direction for all values of T . The 
asymmetric case ( [i = 1.0) in the absence of space dependent friction give current in the positive direction only as 
a function of temperature. Separately in both these cases absolute value of current exhibits a maxima as a function 
of T, reminiscent of stochastic resonance phenomena. In a purely asymmetric case (A = 0) current vanishes rapidly 
when T exceeds the temperature associated with the barrier height. Whereas, in the symmetric case due to space 
dependent friction absolute values of currents are significantly higher and decay slowly to zero in the large temperature 
regime. Naturally in the presence of both asymmetry and space dependent friction for the case under study the low 
temperature regime is dominated by the effect of asymmetry while the high temperature regime is dominated by space 
dependent friction. From this , one can qualitatively explain the current reversals from positive to negative side as 
a function of temperature even in the adiabatic limit. In the regime (f> > n current due to space dependent friction 
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and potential asymmetry flows in same direction and hence no reversal is possible as a function of T. For frequencies 
higher than the interwell frequency wq, the low temperature scenario changes. The direction of current in this regime 
is more of a interplay between potential asymmetry and w than that of A. Due to higher frequency the Brownian 
particles do not get enough time to cross the right barrier which is at a larger distance from the minima. Number of 
particles moving about the potential minima increases. This fact is amply reflected in Fig. |^, where the time averaged 
probability curve P av (x) = i JJ" P(x,t)dt is plotted as function of x for various values of w with T — 0.05. ft is to 
be noted that the distribution is independent of rj(x). Figure ^| shows that the probability of finding the Brownian 
particles near the minima of the potential well increases with increasing frequency, consequently the probable number 
of particles near the potential barrier decreases. Since the distance from a potential minima to the basin of attraction 
of next minima is less from the steeper side than from the slanted side, hence in one period the particles get enough 
time to climb the potential barrier from the steeper side than from the slanted side, resulting in a negative current. On 
increasing the temperature , the particles get kicks of larger intensity and hence they easily cross the slanted barrier, 
resulting in a current reversal and positive current. On further increasing the temperature, the effect of A dominates 
as mentioned previously and as a result the current again flows in the negative direction, implying a second current 
reversal as shown in dotted curve of Fig [j]. It is obvious from the above argument that for higher frequency the first 
reversal temperature will be higher as shown in the Fig. ^. The dotted line in the base of Fig. ^ shows the contour 
of zero current. But as mentioned previously, since the effect of A dominates for higher temperature, therefore the 
second reversal temperature decreases with increasing frequency. Beyond the critical frequency w c there is no reversal 
of current as shown in Fig.( |], §)■ In the high frequency regime the effect of space dependent friction dominates the 
nature of current. We would like to emphasize here that in absence of asymmetry current reversals does not take 
place. 

Multiple current reversals can also be seen when the amplitude ( A) of the forcing term is varied in a suitable 
parameter regime of our system. In Fig. ^, the plot of j versus A is shown for different values of w, keeping A, cf> and 
T fixed at 0.1, 0.887T and 0.05 respectively. For w = 4.0 curve, we can see as many as four current reversals. For very 
large value of A, the current asymptotically goes to a constant value depending on the value of </>, as was previously 
shown for the adiabatic case [0. This value was shown to be — w sin(</>). This is special to the space dependent friction 
ratchet where the currents saturate to a finite value in the large A limit. In the absence of space dependent friction 
it is to be noted that currents decay to zero in the same asymptotic regime. As the asymptotic value depends on the 
phase tf>, so we can choose it appropriately to make the asymptotic current positive or negative. In the present case 
(j) is chosen such that the asymptotic current is negative which guarantees at least one current reversal irrespective of 
frequency. The oscillatory behavior in the j — A characteristics is the reminiscent of the deterministic dynamics [ff|]29|] 
which will be discussed later. The inset in Fig. || shows the zero contour of j versus <f> and A for w = 4.0. For <fi > ir 
only two current reversals can be seen. For it > (f> > ir — e (where e is a small number ) or < <f> < e, four or more 
current reversals can occur. The value of e depends critically on T, large for small T and vice-versa. It is also to be 
noted that we have even number of reversals for finite frequency driving. These oscillatory features along with their 
associated current reversals disappear in the high temperature regime as expected. Figure || is for w — 0.25 and shows 
that current reversal is also possible in the deterministic regime of the overdamped system. This deterministic reversal 
of current cannot be attributed to the chaotic motion of the system like in an forced underdamped oscillator [po[ . 
It solely arises due the presence of space dependent friction. As shown in Ref. here too the deterministic current 
shows quantization and phase locking behaviour. However all these complex features are not robust in the presence 
of noise as discussed in earlier literature. 

In conclusion, we have studied the transport properties of overdamped Brownian particles moving in an asymmetric 
potential with space dependent friction coefficient and rocked by periodic force. We observe several novel and complex 
features arising due to the interplay between asymmetry and inhomogeneous friction. Currents in the low temperature 
regime is mostly influenced by the asymmetry of the potential. At higher temperatures it is controlled by the 
modulation parameter A of the friction coefficient. We find current reversal with temperature even when the forcing is 
adiabatic. In the presence of finite frequency, twice current reversals are seen. As function of amplitude of the forcing 
term we observe multiple current reversals. Current even reverses its sign in the adiabatic deterministic regime. All 
the above results can be understood in a qualitative manner. We expect that our analysis should be applicable for 
the motion of particle in porous media and for molecular motors where space dependent friction can arise due to the 
confinement of particles. 

MCM acknowledges partial financial support and hospitality from the Institute of Physics, Bhubaneswar. MCM 
and AMJ acknowledgepartial financial support from the Board of Research in Nuclear Sciences, DAE, India. 



3 



[1] S. Leibler, Nature 370,(1994),412; J. Maddox, ibid 365(1993), 203; ibid, 368, (1994), 287; ibid, 368, (1994), 287. 
[2] M.O. Magnasco, Phys. Rev. Lett. 71, 1477, (1993); I. Derenyi and T. Vicsek, Phys. Rev. Lett. 75,(1995),374. 
[3] P. Jung, J. G. Kissner and P. Hanggi, Phys. Rev. Lett, 76(1996), 3436. 
[4] R. Bartussek, P. Hanggi, and J.G. Kissner, Europhys. Lett. 28, 459, (1994). 

[5] J. Prost et al, Phys. Rev. Lett. 72, (1994), 2652; J. F. Chauwin, A. Ajdari and J. Prost, Europhys. Lett. 27,(1994), 421; J. 
Rousselet et al, Nature, 370, (1994), 446; C. Van den Broeck et al, Lecture Notes in Physics: "Statistical Mechanics of 
Biocomplexity", Vol527, (Springer- Verlag Berlin, Heidelberg, New York), (1999), Page 93; P. Hanggi and 
R. Bartussek, Nonlinear Physics of Complex Systems - Current Status and Future Trends, Lecture Notes in Physics, 
Vol. 476, ed. by J. Parisi, S.C. Mueller, and W. Zimmermann (Springer, Berlin, 1996), pp. 294-308. 

[6] A. Ajdari and J. Prost, Europhys. Lett. 32, (1995), 373. 

[7] P. Reimann, R. Bartussek, R. Hauss ler and P. Hanggi, Phys. Lett. A, 215, (1996), 26. 
[8] C.R. Doering, W. Horsthemke, and J. Riordan, Phys. Rev. Lett. 72, 2984, (1994). 

[9] A. Ajdari, D. Mukamel, L. Peliti, and J. Prost, J. Phys. (Paris) 14, 1551, (1994); M.C. Mahato, and A.M. 

Jayannavar, Phys. Lett. A209, 21, (1995); D.R. Chialvo, and M.M. Millonas, Phys. Lett. A209, 26, (1995). 
[10] M. M. Millonas, Phys. Rev. Lett 74, 10 (1995). 
[11] A. M. Jayannavar, Phys. Rev. E 53, 2957 (1996). 

[121 M. C. Maha to, T. P. Pareek and A. M. Jayannavar, Int. J. Mod. Phys B 10, 3857 (1996). |cond- 

mat/9603103| . 
[13] M. Biittiker, Z. Phys. B 68, 161 (1987). 

[14] N. G. van Kampen, IBM. J. Res. Develop 32, 107 (1988). 
[15] Rolf H. Luchsinger, Phys. Rev. E, 62, (2000) 272. 

[16] D. Dan, M. C. Mahato and A. M. Jayannavar, Int. J. Mod. Phys. B, (in press); Cond-Mat, (0006244) 
[17] C. Kettner et al, Phys. Rev. E, 61, (2000), 312. 

[18] D. Dan, M. C. Mahato and A. M. Jayannavar, Phys. Lett. A 258, 217 (1999). 
[19] D. Dan, M. C. Mahato and A. M. Jayannavar, Phys. Rev. E 60, 6421, (1999). 
[20] Jose L. Mateos, Phys. Rev. Lett, 84, (2000), 258. 
[21] P. Reimann Phys. Rep., 290, (1997), 149. 

[22] M. Bier and R. D. Astumian Phys. Rev. Lett., 32, (1995), 373. 

[23] Luc P. Faucheux and A. J. Libchaberm, Phy. Rev. E, 49, (1994), 5158. 

[24] H. Brenner, Chem. Eng. Sc., 16, (1962), 242. 

[25] C. M. Falco, Am. J. Phys., 44, (1976), 733. 

[26] A. M. Jayannavar and M. C. Mahato, Pramana- J. Phys 45, 369 (1995). 

[27] J.M. Sancho, M. San Miguel, and D. Duerr, J. Stat. Phys. 28, 291, (1982). 

[28] H. Risken, The Fokker Planck Equation (Springer Verlag, Berlin , 1984). 

[29] T. E. Dialynas, K. Lindenberg and G. P. Tsironis, Phys. Rev. E, 56, (1997), 3976. 



4 



0.0025 



0.0015 



0.0005 

■ 

J 

-0.0005 



-0.0015 



-0.0025 



co = 3.0 
co = 4.0 
oo = 5.0 





0.0 0.1 0.2 0.3 0.4 



0.0 



0.2 



T 



0.4 



0.6 



FIG. 1. The mean current j vs temperature T for cj> = 0.27T, A = 0.5 and A = 0.1. The driving frequencies are w = 3.0, 4.0 
and 5.0. The inset shows the contour of zero current for same values of T, A. Regions enclosed on the right hand side of a given 
contour is the negative current region and vice versa. Note that we have only one reversal for all values of cf> > it. 
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FIG. 3. The time averaged probability P av vs x for various values of w and T = 0.05 and A = 0.5. The distribution is 
independent of <f> and A. 
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FIG. 5. The mean current j with amplitude A of the forcing term for <j) — 0.887T, T = 0.05 and A = 0.1 with w = 3.0, 4.0, 5.0. 
The inset shows the contour of zero current for T = 0.05 for w = 4.0. The dots in the inset shows the four values of A for 
4> = 0.9n where the current reversal occurs. 
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FIG. 6. The deterministic current j vs amplitude of forcing A for w = 0.25 and <j> = 0.2n and A = 0.1. 
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